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ABSTRACT 


The responses of an unsymmetrically laminated fiber-reinforced composite 
cylinder to an axial compressive load, a torsional load, and the temperature 
change associated with cooling from the processing temperature to the service 
temperature are investigated. These problems are consider axisymmetric and 
the response is studied in the context of linear elastic material behavior and 
geometrically linear kinematics. Four different laminates are studied: a gen- 
eral unsymmetric laminate: two unsymmetric but more conventional laminates: 
and a conventional quasi-isotropic symmetric laminate. The responses based 
on closed-form solutions for different boundary conditions are computed and 
studied in detail. Particular emphasis is directed at understanding the influ- 
ence of elastic couplings in the laminates. The influence of coupling de- 
creased from a large effect in the general unsymmetric laminate, to practically 
no effect in the quasi-isotropic laminate. For example, the torsional loading 
of the v general unsymmetric laminate resulted in a radial displacement. The 
temperature change also caused a significant radial displacement to occur 
near the ends of the cylinder. On the otherhand, the more conventional un- 
symmetric laminate and the quasi-isotropic cylinder did not deform radially 
when subjected to a torsional load. From the results obtained, it is clear the 
degree of elastic coupling can be controlled and indeed designed into a cyl- 
inder, the degree and character of the coupling being dictated by the applica- 
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Chapter 1 - Introduction 


Since their introduction as a structural material, 
fiber-reinforced composites have been used primarily as a 
replacement for more traditional materials such as metal. 
Despite the wide range of tailoring, elastic couplings, and 
directional dependence available with fiber-reinforced 
composite materials, they have not been used with these 
characteristics in mind. Rather, the inherent lighter 
weight, integrated design, and near-net-shape manufacturing 
capability of these materials have been the prime reasons 
for their increased use. As a result, composite structures 
today are overweight, conservatively designed, and in many 
applications, not that much better than, say, a lithium 
aluminum structure designed for the same purpose. There 
are exceptions to this, the forward sweep experimental 
aircraft technology being a prime example. While in some 
applications a design that does not utilize fiber- 
reinforced materials may be the best design, it is fair to 
say at this point in the development of composite 
structures that their potential in many applications, 
traditional and not so traditional, has not been fully 
explored. It has not been fully explored for two major 
reasons. First, the majority of thinking is still metals, 
or isotropic, in nature. A culture based on the 

availability of directional dependency and elastic 



couplings does not yet exist. Second, analysis of 
structures and materials with nonisotropic behavior and the 
wide variety of elastic couplings available with composite 
materials has not matured. For sure it is possible to use 
numerical methods such as finite-elements and obtain 
numbers for a particular problem, or a series of problems. 
However, such an approach provides no insight as to what 
characteristics of composites materials are important and 
may provide benefits to a particular class of problems. 
Solving simple problems in the context of the more complex 
behavior of composite materials can lead to insight into 
how certain characteristics of the responses of composite 
materials influence the behavior of structures. This 
latter approach to understanding the behavior, and hence 
understanding the full potential of fiber-reinforced 
composite materials, is one of the purposes of this thesis. 
Specifically, this thesis examines the response of 
cylinders fabricated with fiber-reinforced materials, the 
cylinders being presumably filament wound from a 
graphite/polymer, and the winding sequence being completely 
general, in particular, unsymmetric. Herein the term 
'unsymmetric' will be used as the commonly accepted 
description of a composite laminate that does not exhibit 
elastic symmetry with respect to its geometric midplane. 
It is well known that due to the elevated temperatures 
associated with the processing of polymer matrix 


composites, an unsymmetric laminating sequence in a flat 
laminate can cause severe warpage of the laminate as it 
cools to a lower service temperature [1, 2]. An 
unsymmetrically laminated cylinder, on the other hand, 
will remain cylindrical when cooled from the elevated 
processing temperature. Thus, unsymmetrically laminated 
cylinders offer a real possibility of exploiting the 
advantages of unsymmetric laminates, if there are indeed 
any advantages. Clearly what may be considered an 
advantage depends on the application. Hence, understanding 
and cataloging the response of unsymmetrically laminated 
cylinders is an important first step, before specific 
applications are investigated. Some work has been done in 
this regard but only for one very special case. [3] 

Here the response of an unsymmetrically laminated 
cylinder to an axial compressive load, a torsional load, 
and due to the temperature change associated with cooling 
from the processing temperature to the service temperature 
are investigated. These problems are all considered 
axisymmetric, and the response will be studied in the 
context of linear elastic material behavior and 
geometrically linear kinematics. Clearly these studies are 
the precursor to studies aimed at investigating the 
potential for utilizing unsymmetrically laminated 
construction for delaying the collapse of cylinders due to 
either an axial load, a torsional load, or both. Such 


problems are geometrically nonlinear problems and can only 
be studied after completely understanding the linear 
counterparts. The thermal problem is important because 
there may be thermally induced warpage of the cylinder near 
its ends. Such warpage would influence the pre-collapse 
response in the nonlinear analysis. The second purpose of 
the work described here, then, is to develop and solve the 
equations which govern the behavior of unsymmetrically 
laminated cylinders subjected to axial compression, 
torsion, and a temperature change, and translate the 
solutions into software so numerical results for specific 
examples can be easily obtained and studied. 

The next chapter introduces the basic nomenclature 
relevant to fiber-reinforced composites as it pertains to 
the present class of problems. Emphasis is on the 
constitutive behavior of the material. The third chapter 
develops the equations assumed to govern the behavior of 
the cylinders, namely, Sanders' kinematic approximations 
[4] in conjunction with the constitutive behavior and the 
first variation in the total potential energy. [5] The 
variation formulation is preferred over the Newtonian 
approach because boundary conditions associated with the 
kinematics of the problem are developed as part of the 
variational process. The proper boundary conditions are 
not as easily of ormulated when using the Newtonian approach. 

The fourth chapter briefly outlines the computer 


programs written to implement the closed-form solution of 
the governing equations for the various cases studied. The 
fifth chapter presents numerical results for the three 
loading cases, namely compression, torsion, and a 
temperature change. The responses of several laminates are 
studied, one laminate being more unorthodox than the 
others. The unorthodox laminate is studied to explore the 
full range of responses that can be exhibited with 
composite materials. Two laminates felt to be viable 
alternatives to current conventional designs are also 
studied. These laminates are compared to what could be 
considered a conventional counterpart. Concluding remarks 
and recommendations for future research are discussed in 


Chapter 6. 



Chapter 2 - Some Basic Definitions Associated with 


Fiber-Reinforced Composite Materials 


A single layer of a fiber reinforced composite 

o 

material is referred to as a lamina. The lamina coordinate 
system, illustrated in fig. 1, defines the x axis in the 
fiber direction, the y axis perpendicular to the fibers and 
in the plane of the lamina, and the z axis perpendicular to 
the plane of the lamina. For an anisotropic material 
Hooke's Law can be defined as 
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This results in 21 independent elastic constants, c ij ' 
being required to describe the material. The C matrix is 
called the stiffness matrix. By considering materials 
which exhibit some degree of symmetry, the number of 
independent elastic constants is reduced. For example, 
particular to the study here, an orthotropic material has 
only 9 independent constants. Hooke's law then takes the 
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The stiffness matrix can be written in terms of the 
engineering material properties, namely Young's moduli, 
Poisson's ratios, and shear moduli. 

For the plane stress state, a z = t xz = Ty Z = 0, and 
the 6x6 matrix can be reduced to a 3 x 3 matrix. For this 
situation Hooke's Law for an orthotropic material becomes 
[ 6 ] 
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The Q matrix is called the reduced stiffness matrix and is 


defined as 
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where E x and Ey are Young's moduli in the x and y 
directions, respectively; n x and p,y are Poisson's ratios 
when the material is loaded in the x and y directions, 
respectively; and E s is the shear modulus of the material. 
The reduced stiffness matrix can obviously be related to 
the Ci j . Here it is more convenient to use the engineering 
properties directly. 

When several laminae are bonded together they form a 
laminate. Each lamina of the laminate can have a unique 


9 



fiber angle relative to a global, structural, or laminate 
coordinate system. For example, a cylinder may have one q 

layer wound helically at 45* to the cylinder axis, another 

layer wound at -45* to the cylinder axis, and yet a third 
layer with fibers parallel with the axis. Since any # 

analysis is usually conducted in the global system, the 
reduced stiffness matrix for each lamina must be 

transformed from the lamina coordinate system to the global 
or laminate coordinate system. The global coordinate 
system is denoted as the 1-2 system and the transformation 
is shown in fig. 2. This transformation results in what is 
called the transformed reduced stiffness matrix. Thus, for 
the plane stress case, referred to the laminate coordinate 
system, Hooke's Law becomes 



with the Q matrix defined as 



( 10 ) 


( 11 ) 
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j where m = cos 0 and n = sin 9, 0 being defined in fig. 2. 

The subscripts 1 and 2 on the stresses and strains in eq. 
10 indicate these quantities are referred to the 1-2 
coordinate system. 

i 
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Assuming the Kirchhoff hypothesis is valid for a 
laminate, [5,6] the strains, €]_, e 2 , e 12 , of eq. 10 can be 
written in terms of the laminates's midplane extensional 
strains, e*i, e* 2 , e *12' an< * the midplane curvatures, k 0 ^, 
k *2' k *12' as 
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The stress resultants for a laminate can be defined as 
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pt/2 

M. _ a J z t. . dz , (25) 

12 z=-t/2 12 

where t is the thickness of the laminate and it is implied 
that the z = 0 location is at the laminate geometric 

midplane. The N's are the force resultants and the M's are 
the moment resultants. 

Since the midplane strains and curvatures are not 
functions of z, substitution of eq. 10 into eq. 20 results 
in 



Because fiber angles can vary from layer to layer, the Q's 
do vary with z but are piecewise constant with z. The 
integrals become sums, the sums being defined as 


14 




(28) 
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where z^ is the z location of the interface between the k*-* 1 
and the (k-l) t * 1 layers, and z 0 * -t/2 and z n = +t/2. 

Similar relations can be derived from eqs. 22 and 24 
to produce a 3 x 3 A matrix. The A matrix is called the 
laminate inplane stiffness matrix and an abbreviated form 
can be written as 
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Considering egs. 20, 22, and 24, the B^j defines another 
3x3 matrix. The B matrix is called the laminate 
bending/stretching stiffness matrix. 

The results so far lead to a relation of the form 
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More will be said regarding the B matrix shortly. 

The moment equations, eqs. 21, 23, 25, can be defined 
similarly. For example, substituting eq. 10 into eq. 21 
results in 
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since the Q's are piecewise constant with z, the integral 


17 



becomes a sum, namely 



N 


layer 





Considering also eqs. 23 and 25, the Djj forms a third 
3x3 matrix. This matrix is referred to as the laminate 
bending stiffness matrix. 

Arrangement of the six stress resultants into matrix 
form leads to the load-midplane deformation relationship 



(36) 


The 6x6 matrix of eq. 36 is called the laminate 
stiffness matrix, or the ABD matrix. When viewed in the 
context of eq. 36, it is clear that composites exhibit 
elastic couplings simply not present with isotropic 



L8 


materials. The existence of the elements of the B matrix 
implies that force resultants (N^ , N 2 , N 12 ) and curvatures 
are coupled, and moment resultants (M lr M 2 , M 12 ) and 
extensioned strains are coupled. Because of the D 16 and 
D 2 6 terms, bending moment resultants (M^, M 2 ) cause twist 
curvatures (k^*) and twisting moment resultants (M 12 ) 
cause bending curvatures (k^*, k 2 * ) . This is a very 

powerful characteristic of fiber-reinforced materials. For 
specific fiber orientations, some of the elements of the 
ABD matrix are zero. For a symmetric laminate B is 
identically zero. Furthermore, if for every layer with 
fiber angle 6 there is an identical layer with fiber angle 
-6, the laminate is said to be balanced and A^g and A 26 are 
zero. (See fig. 2.) Finally, if each layer has a fiber 
angle of 0* or 90*, D^g » D 25 = 0. While the analysis of 
composite structures is simpler if these elements of the 
ABD matrix are zero, it is not clear there is any 
structural advantage to having these zero. It is the 
specific intent of this thesis to investigate the 
ramifications of them not being zero. 

Writing out eg. 36, 
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Thermal Expansion Effects 


Thermal expansion effects of the fiber-reinforced 


20 


material can be included by returning to Hooke's Law. With 
thermal expansion effects Hooke's Law becomes 



(43) 


where a^, a 2 , and a 32 are lamina coefficients of 
thermal expansion as measured in the global 1-2 coordinate 
system. The quantities a^, a 2 , and a 3 are defined as 
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where a x is the lamina coefficient of thermal expansion in 
the fiber direction, and ay is the lamina coefficient of 
thermal expansion perpendicular to the fibers in the plane 
of the lamina. 

Using the definitions of force and moment resultants, 
eqs. 20-25, and substituting for a 2 , and t 12 from eq. 
43, each force and moment resultant can be calculated to 


(44) 
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include thermal expansion effects. The integrals which led 
to the ABD matrix are not affected by the inclusion of the 
thermal effects. However, the thermal effects do lead to 
additional integrals, with thermal effects, for example, 
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Since it has the proper units, the integral is interpreted 
as an equivalent thermal inplane load, N^, i.e., 
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Since the Q's and a's are piecewise constant with z, the 
integrals can be replaced by summations, namely, 
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Considering N 2 , N i2' M l' M 12 results in similar 

definitions of thermally-induced loads, namely. 
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( 53 ) 


M 


12 


N layer 

U ‘X % + % + Q66 - ’ 


’k X *k 


* - 7 . * 


( 2S 

1 2 


'k-i 


) AT . 


(54) 


The latter three are thermally-induced moments. 

With the above developments, the final form from the 
stress resultants, including the thermal effects, is 


N 1 * A 11 e i + A 12 e 2 + A 16 6 12 


+ B X1 k* + B 12 k* + b 16 k£ 2 


- N 


N 2 ~ A 12 6 i + A 22 £ 2 + A 26 e l2 


(55) 


+ B 12 k l + B 22 k 2 + B 26 k 12 " N 


2 

( 56 ) 


N 12 ~ A 16 € 1 + A 26 6 2 + A 66 £ 12 
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* 


* 


+ B 16 k i + B 26 k 2 + B 66 k 12 


- N 


12 

( 57 ) 


M. 


B 11 6 1 + B 12 € 2 + B 16 £ 12 


+ k* + D 12 k* + D lg k£ 2 


- M 


1 

(58) 


M, 


B 12 e i + B 22 e 2 + B 26 e l2 


°12 k i + °22 k 2 + °26 k 12 


- M 


2 

(59) 


M 12 = 


B 16 e i + B 26 e 2 + B 66 e i2 


+ °16 4 


k° + k* - M, _ . 


26 2 "66 12 


12 
(60) 


The definition of the ABD matrix and the inclusion of 
thermal effects are important ingredients for the next 
step, namely, the development of a structural level theory 
for cylinder behavior. A quick review of this chapter 
indicates that all that has been done is to define 
thickness-integrated properties. These relations developed 
shift emphasis to the laminate midplane, yet nothing has 
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4 

been said as to how the N's, M's, e's, or k's vary from 
point to point on the midplane, say, with distance along 
the cylinder. This is done in the next chapter. * 
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Chapter 3 - Sanders' Shell Theory: 
Axisvmmetric Problem 

3 . 1 The Variational Statement 

To derive the equations governing the deformation of 
a cylindrical shell, the variational method in combination 
with the total potential energy is used. The total 
potential energy of a three dimensional solid is [5] 


r * U - J J J v B. U. dV - J T{ V) u i ds, 

(61) 

where U is the strain energy of the solid and the other 
terms are due to the applied loads. Since this study will 
not involve applied loads, but rather displacement boundary 
conditions, 


J"! I v B. Ui dV - J J s T‘ V) Ui 


dS = 0 (62) 


and so 


r = u 


(63) 


With this, the first variation of the total potential 
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energy reduces to 



: °x ie x + 


°e Se e 


S e. 


T er 5e er + 


t Se + 
xr xr 


T xe 5e xe > dv - 


(64) 


This expression will be applied to cylinders. Figure 3 
illustrates the coordinate system used to describe a 
cylinder with a mean radius R and a length L. The cylinder 
thickness will be denoted as t, and, as in the last 
chapter, t represents the thickness of the laminate which 
constitutes the cylinder wall. The subscripts x, 0, and r 
on the stresses and strains in eq. 64 indicate that these 
quantities are referred to the global cylindrical 
coordinate system. 


3 . 2 Kinematic Assumptions 

To use the variational statement of eq. 64, the 
strain field must be defined. Since the strains are 
related to the displacements, it is the displacement fields 
that must be defined. For the axisymmetric case, the axial 
displacement, u(x,z), and the radial displacement, w(x,z), 
of an arbitrary point A are illustrated in fig. 4 and are 



FIGURE 3 THE CYLINDER COORDINATE 
SYSTEM 


6Z 



X 



DEFORMED 

CROSS-SECTION 



FIGURE 4 THE AXIAL AND RADIAL 
DISPLACEMENTS OF 
CYLINDER 
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equal to 


(65) 

( 66 ) 

where the superscript * refers to the displacements of the 
cylinder reference surface which is located at the mean 
radius. The distance z is measured from the reference 
surface, i.e., z = 0 is the geometric midsurface of the 
cylinder wall. The circumferential displacement, v(x,z), 
of an arbitrary point A is illustrated in fig. 5 and it is 
found to be 

v(x,z) = (1 + | ) V (x) . (67) 

This form of v, as opposed to 

v ( x , z ) = V* (X) , (68) 

is referred to as the Sander's kinematic assumption [4], 
while the latter is attributed to Donnell [7]. With 

Sander's kinematic assumption, the validity of the 
Kirchhoff hypothesis is implied. Note: The three problems 

of compression, torsion, and temperature change are assumed 
to be axisymmetric and so the displacements are not 


u(x, z) = u # (x) - z 


dw’ (x) 
dx 


W(X,Z) * W* (X) , 
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DEFORMED 



FIGURE 5 KINEMATICS OF CIRCUMFERENTIAL 
DISPLACEMENTS OF CYLINDER 
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functions of 6. Furthermore, the stresses, strains, and 
the stress resultants are assumed to be independent of 0 
and so derivatives with respect to 0 will be set equal to 
zero. 

Knowing the displacements, the strain field can now 
be defined using the strain-displacement relations. 


3.3 The Strai n Field 


The strain-displacement relations in polar 
coordinates for the axisymmetric case are as follows [8]: 


e x (x,r) = 


^u(x.r) 

^ x 


(69) 


e 9 (x, r) = 


w(x,r) 

r 


(70) 


e r (x,r) - 


^w(x, r) 
<*r 


(71) 


6 er (x ' r) 


- v ( x , r ) 
r 


<^v(x.r) 

or 


(72) 


e xr (x ' r) 


^u(x,r) 

ox 


+ <^w(x, r) 

ox 


(73) 
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♦ 


e xe (x ' r) 


ax 


(74) 


Substituting the displacements of eqs. 65-67 into the 
strain-displacement equations, eqs. 69-74, leads to 


« x (x,z) 


du* (x) 
dx 


d t w* (x) 
dx* 


(75) 


e 0 (x,z) 


W (x) 
R 


(76) 


e r (x,z) 


d W(x) = 0 

dz 


(77) 


e 9r (x ' z) = 0 


(78) 


e xr (x ' 2) = 0 


(79) 


6 xe (x ' z) 


dv‘ (x) z dv* (x) 

dx R dx 


(80) 


♦ 


4 


where in the above equations it has been assumed that the 
shell is thin so 


r = R+ z = R (1 + |)*R 


(81) 


and use has been made of the fact that 


4 
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In addition, the functional dependence on r has been 
switched to the functional dependence on z. As can be 
seen, three of the six strains are zero, the nonzero ones 
, being e x (x,z), e 0 (x, z) , and e x0 (x,z). Therefore, the 
variational statement, eq. 64, simplifies to 

5r - Jv < °x U X + °6 S£ 8 + T x© S£ xe > dV - < 83 > 


3 . 4 The Euler Equations and Boundary Conditions 


Substitution of 

eqs. 75, 

76, 

and 

80 into 

the 

simplified variational 

statement, 

eq. 

83, 

leads to 

the 


governing equations and boundary conditions. The quantity 
dV in polar coordinates is approximated for thin shells as 
follows: 


dV = 

r d© 

dr 

dx 



(84) 

dV ■ 

(R + 

2) 

de 

dz 

dx 

(85) 

dV = 

R ( 1 -t 

■ z/R) 

de 

dz dx * R d© dz dx. 

(86) 


Using eqs. 84-86 and eqs. 75-80, the variational statement 



of eg. 83 becomes 


5 r = 


jL/2 

X=-L/2 


I 

8—0 


2ir 


J V2 

Z=-t/2 


t a x « ( 


du* fx) 
dx 


, d^w* _ ,, w* , , dv* 

z J..1 ) + <?Q o ( p ) + T vfl ° ' rtv 


dx ! 


x8 v dx 


+ R-S‘> ] dz Rde dx ' 


(87) 


where the limits on the integral define the volume of the # 

cylinder, and the points x =* 0 and z = 0 are located at the 
mean radius midway along the cylinder length. Integration 
with respect to z can be done independently and in doing so % 

the following stress resultants are defined: 


p t/2 

N x 3 J °x dz 

X Z»-t/2 X 


pt/2 

M a J z a dz 

Z=-t/2 


pt/2 

N 0 » J a dz 

0 Z=-t/2 0 


( 88 ) 


(89) 


( 90 ) 
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< 


(91) 


N 


I 


t/2 


vQ — j T V0 *^ Z 

X9 Z=-t/2 X0 


M 


I 


t/2 


ye — J Z T V0 • 

xe Z=-t/2 xe 


(92) 


With the integration on 0 and substitution of eqs. 88-92, 
the variational statement of eq. 87 becomes 


sr 


r l/2 

2itR J [ N 

x=-L/2 3 


c, 

5( dx } 


H 2 w * 

"x 8 < ZF I 



5 “‘ + N xe 5 < S‘> 




+ 


M x0 c , dv‘ 
R 5 ( dx 


dx. 


(93) 




Integration by parts leads to 


5r = 2-jrR 


N x SU- 


L/2 

-L/2 


( M x S 



L/2 

-L/2 
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+ ( -di 5w- >i 


-L/2 


+ (N xe 5v ’> 


-L/2 


M X0 V 2 

+ ( -§2 SV) | 

R -L/2 


pL/2 ”dN _VA11 xe 

1 t ( ) Su ’ + ( -d? 


X=-L/2 


M 

d f xQ. i c„< 

dx ( R > }5v 


-d*M v N 

+ { + “f > *W ] dX 


(94) 


Of course the first variation must be zero and so eq. 94 
reveals the Euler equations and the associated boundary 
conditions for the problem. The Euler equations are: 


(95) 


& < N xe + 


( 97 ) 


while the boundary conditions are: 
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0 or 


u° specified 


(98) 


either N x = 


, (}m • 

either M x = 0 or ^ specified 


(99) 


either 


dMj 

dx 


= 0 or w # specified 


( 100 ) 


M. 


X0 


either ( N x0 + -jf ) 


0 or v* specified. (101) 


It should be noted that the quantities N x0 and M x0 /R always 
appear summed. This quantity will be referred to as N x0 , 

X • 6 • f 



N 


X0 



(102) 


and is interpreted to mean that the Sanders ' theory does 
not use either N x0 or M x0 separately. Donnell's theory, on 
the other hand, makes no mention of M x0 . Obviously M x0 /R 
is a result of the z/R term in eq. 67. 


3.5 Solving the Euler Equations 


Referring to eqs. 17-19 and 75-80, and using the 
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notation of Ch. 2 by letting the subscripts 1 = x and 2 = 
0, for the cylinder problem 


e x 


du‘ 

dx 




w* 

R 


€ xe 


dv * 

dx 


d^w' 

dx* 


v • 
*0 


k • 
K x0 


1 dv ‘ 
R dx 


(103) 


(104) 


(105) 


(106) 


(107) 


(108) 


Rewriting eq. 36 for the x-0 coordinate frame, substituting 
eqs. 103-108 into eq. 36, and writing out each stress 
resultant with the equivalent thermal loads included leads 
to the expressions for the stress resultants pertinent to 
axisymmetric cylinder problems. These stress resultants 
are: 


N 


x 





w° 

R 




dv" • 

dx 
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(109) 




- B 


dV 
11 dx* 


T 

N x 


N = A + A — * + 

w 0 12 dx *22 R 


< A 26 + R B 26 ) £' 


d*w * 
12 dx* 


- B 


N 


6 


( 110 ) 


N 


X0 


(A + hi > du' 

1 *16 R } dx 


• / a 4 . ^26 w* . , 2 B 66 . °66 . dv‘ 

( A 26 + R ) R + ( + 13 + o* ) 


‘66 


R* ; dx 


/ B + 116 dlw* 

( B 16 + R } dx* 


N xe - < lll > 


M. 


B du* + B w* , 1 n x dv' 

B 11 dx B 12 R B 16 R °16 dx 


d*W 

11 dx* 


- D 


M 


( 112 ) 


M e - B 12 + B 22 I' + < B 26 + k °26 » 


- D 


d*W 
12 dx* 


M 


0 


(113) 


where 
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(114) 


Using eqs. 109-114, the Euler equations, eqs. 95-97 


A + (A 

A 11 dx* ( 16 


B x 1 . dw‘ 

°16 ' dx* R 12 dx 


d 3 w * 
11 dx 3 


(115) 


. 1 R \ dju* 

( A 16 + R B 16 * dx* 


< A 66 + i B 66 


I n ) 4!v' 

R* °66 ' dx* 


(A + — B ) — 

{A 26 R °26 ; R dx 


, B + I D x dV 

‘ 16 R 16 ’ dx 3 


(116) 


11 dx 3 


I A du' + , B + I D \ d^v‘ 

R A 12 dx + B 16 R °16 ' dx 3 


42 


R < A 26 + 


I B 

R B 26 


\ dv’ 

' dx 


- D 


d 4 w ° 
11 dx 4 


1 B 

R B 12 


d*w * 

dx* 


R- A 22 “• 


1 T 

— M 

R 0 


0 


(117) 


Before proceeding, some of the influences of elastic 
couplings should be discussed in the context of the stress 
% resultants, eqs. 109-113, and the governing equations, eqs. 

115-117. For a symmetric laminate the B^j are equal to 
zero and the stress resultants equations become 


N = A — + 

x A 11 dx 


12 


w* 

l + 


. dv* 
A 16 dx 



(109a) 


N e - A i2 S' + 


W * 

A — + 

A 22 R 


dv* 
26 dx 



(110a) 




du* 

dx 


+ 


A - 
26 R 
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i 


( a, + -££ ^ £lv* 

+ ^ A 66 + R* ^ dx 


16 d *w ‘ 
R dx* 


T 

N xe ' 


M. 


I D 32 * 

R U 16 dx 


d*w ‘ 
’ll dx* 


M 


e 


I D 

R u 26 dx 


. d*w ' 
i2 dx* 


The differential equations become 


, dlu* , , dW* 1 . dw' 

A 11 dx* 16 dx* R a 12 dx 


0 , 


A + 

A 16 dx* 


( A 66 + °* ^ 


I n ) Hi* 

R* °66 ) dx* 


1 dw‘ 
^26 R dx 


I D 

R 16 dx 3 


o , 


i 

(Ilia) 

I 

(112a) 

i 

(113a) 

1 


(115a) 


(116a) 
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1 du* 

R A 12 dx 


+ 


R °16 


d 3 v * 

dx 3 


1 dv* 

R A 26 dx 


4 

- D 


d w * 
11 dx 4 



■ 0 , (117a) 

Comparing eqs. 109-113 with eqs. 109a-113a, it is seen that 
with an unsymmetric laminate the curvature of the radial 
displacement in the axial direction, i.e., d'wVdx*, 
contributes to the axial and circumferential force 
resultants through B X1 and B 12 , respectively. For 
symmetric laminates the radial displacement influences the 
axial and the circumferential force resultants only through 
the A ^2 and A 22 , respectively. These A terms are related 
to Poisson effects. Similarly, B i;l and B 12 cause the axial 
displacement gradient, du*/dx, to influence the axial 
moment resultant, M x , and the circumferential moment 
resultant, Mg, respectively. For symmetric laminates the 
axial displacement gradient has absolutely no influence on 
the moment resultants. It is interesting to note that even 
for a symmetric laminate, the circumferential displacement 
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gradient, dV/dx, influences the axial and the 
circumferential force resultants through A 16 and A 2 6 t 

respectively, the shear force resultants through A 66 and 
Dgg respectively, and the moment resultants through D 16 and 

d 26 • t 

For a balanced symmetric laminate, A 16 = A 26 = 0 anc * 
the stress resultants become 


N 


x 


k ll 


du* 

dx 



w* 

R 



(109b) 


i 


XT X du * , 

N 0 A 12 dx 


w* 


22 



(110b) 


N 


xe 



dv* 

dx 


°16 d*W 
R dx 1 


(111b) 


M. 



dv* 

dx 


'll 


d*W 

dx* 


(112b) 


M, 


9 


R °26 


dv* 

dx 



d*w * 

dx* 


(113b) 


While the governing equations become 


I 


46 


i 



(115b) 



d*u “ 

dx* 


+ 


1 dw° 

R A 12 dx 


= 0 , 




d*v * 

dx* 


1 

R 




0 


(116b) 


1 a du* 
R A 12 dx 


+ 


R D 16 


d 3 v * 

dx 3 



. d 4 w * 
11 dx 4 



0 


(117b) 


The axial and circumferential force resultants do not 
depend on the circumferential displacement gradient, 
dV/d>C/ and the shear force resultant does not depend on 
the axial displacement gradient, du*/ d x* For an 
orthotropic laminate and metal, D 16 = °26 = and there is 
further decoupling of displacements and stress resultants. 

In general, the meaning of these couplings is quite 
profound. However, for the cylinder the presence of the 
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couplings just discussed mean that in an unbalanced and/or 
an unsymmetric laminate, twisting causes axial forces (see 
eq. 109) and radial growth or bulging (see eq. Ill) , and 
axial compression causes twisting (see eq. 109) and radial 
growth beyond the normal Poisson effect (see eq. 109) . 


3 . 6 Solution of the Governing Equations ^ 

Equations 115-117 constitute a system of ordinary 
linear differential equations with constant coefficients. ^ 

There is a homogeneous and a particular solution. The 
solution to the homogeneous portion of these equations is 
assumed to be of the form: f 


u* (x) = U e s x ; V(X) = V e s x 

w* (x) * W e s x . (118) 

Substituting these assumed forms into eqs. 115-117 
results in 


[ A 


11 


s 2 ] U 


C (A 


16 


R B 16 ) 


] v 
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+ C I A 12 3 - B 11 s3 J W - 0 ' < 119 > 

t < A 16 + R B 16 ) s 2 1 u 


+ t < A 6S + I B 66 + R* °66 ) s2 1 V 


+ t < A 26 + 


R B 26 > R s - < B 


16 


R D 16 > s 3 1 w 


0 , 


( 120 ) 


t R A 12 3 


B 11 3 I 0 


+ t < A 26 + R ®26 ) R 3 


< B 16 + 5 °!6 > ^ 1 V 


[ D 1X + 5, A 


1 

R* "22 


I B 12 s2 ] W 


( 121 ) 


Writing eqs. 119-121 in matrix form and introducing 
the matrix T^j , leads to 
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Equation 122 is a standard eigenvalue problem and the 
eigenvalues, s^, are found by setting the determinant of 
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the T matrix equal to zero and solving for the s's 
determinant is of the form 


The 


C s 8 + C s 6 + C s 4 = 0 . (129) 


The quantities, Cg, Cg, and C 4 , of eq. 129 are constants 


and are equal to 


2 A U B 66 -t 


2 A 16 B 16 


2 A n b i 6 -r + 


2 A B -ifi 
2 A 16 “ll R 


2 B 


11 R 


+ 2 B 


11 R 


A D 111 

A 11 06 R* 


16 R* 


’ll R* 


11 R* 


2 B 11 B 16 


+ A 11 A 66 D 11 
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( 130 ) 


A 16 D 11 


A 11 B 16 + 2 A 16 B 11 B 16 


A 66 B 11 


2 A 


16 R 


2 A 12 B 16 -Jf 


2 A A -H 
* A 12 ©6 R 


2 A A — =-= 
A 16 a 26 R 


? A A — ix 

Z 11 26 R* 


2 A A — ±-X 
* A 12 A 16 R* 


B 26 B 16 
2 B 11 B 16 “Jf + 4 A 16 B 12 Hr 


12 R 3 


2 A 


12 R 2 


2 A 26 B 11 S? 


2 A 11 B 12 


+ 2 A 12 B 11 "^3 


2 A 11 B 26 
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4 


4 


A,, B 


B 


66 


11 12 R 2 


B 


A, _ B 


66 


12 11 R 2 


+ 2 A, , B 


B 


26 


11 16 R 2 


B 


2 A, - B 


26 


16 11 R 2 


B 


+ 2 A A 

A 11 26 R 


16 


B 


2 A,, A 


16 


12 16 R 


B 


2 A 11 *66 R 


12 


C = A A 

4 A 11 A 22 R 2 


2 ^66 
l 12 R 2 


L 226 

11 R 2 


2 A 12 A 16 


2 6 
R 2 


16 


22 


2 A 


B 


66 


12 


+ 2 A A — — 

2 A 11 a 22 r 3 


2 ^66 
12 R 4 


+ A A 

11 22 r 4 


2 A A ^26 
11 26 _3 


(131) 




Solving eq. 129 for s results in eight eigenvalues. Four 
eigenvalues are repeated zeros and the other four are given 
by 



(133) 

The associated eigenvectors are defined by $ and <p. Here 
U L = *. W. (134) 

V i = 0. (135) 

i = 1,2 ,3, 4, corresponding to the four nonzero roots, s-^. 
To find the proper form for the repeated zero root 
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solutions, assume the following displacements: 


u* (x) = U 5 + U 6 x + U ? x 2 + U 8 x 3 (136) 

V* (x) - V 5 + V 6 X + V ? x 2 + v 8 x 3 (137) 

V (X) = W 5 + W 6 x + W ? x 2 + Wg x 3 . (138) 

Equations 136-138 and their derivatives are substituted 
into equations 115-117. After matching terms with like 
powers of x, it is found that some of the constants in eqs. 
136-138 are zero and what remains is 


u*(x) 

= u 5 + 

U 6 X 


(139) 

v (x) 

+ 

in 

> 

II 

V 6 * 


(140) 

V (X) 

II 

S 

Ul 

• 



(141) 

Substitution 

of eqs. 

139-141 

into eq. 117 

reveals 

relationship 

between W 5 , 

Ug , and 

V 6 such that 
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K R 

w 5 - ( - -qr } u * + 


(-A 


26 


B 26> R 


22 


6 ' 


( 142 ) 


or 


W 5 - a 0 6 + M, 


(143) 


where 


A 12 R 

<- — ) 
A 22 


(144) 


(-A 


26 


B 26> R 


22 


(145) 


By inspection, the particular solution to the system 
of differential equations is 



(146) 


€ 



(147) 




with the subscript p denoting the particular solution. 
Collecting the various portions of the solutions leads to 
the following expressions for the displacements: 
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« 


s l x 


U* (X) = w ! e + 0 2 w 2 e 


S 2 X 


s.x s.x 

+ 0 3 W 3 e + 0 4 W 4 e 


+ U 5 + U 6 x ; 


( 148 ) 


S.X S.X 

V (x) = W, e A + *. W. e 


1 ”1 


2 2 


S 3 X 

* 3 W 3 e 3 + 


*4 W 4 e 


S 4 X 


+ V. + 


V x 
6 x 


(149) 


s.x s.x s.x 

V (x) * W 1 ® + W 2 e + w 3 e 


W 4 e 


S 4 X 


a U 6 + 3 v 6 


-B_ n t 
a 22 N e 


(150) 


As it turns out, the eigenvalues, s^, and assorted 
eigenvectors, and <p in eqs. 148-150 are complex. This 
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complicates the algebra. Since the displacements are 
real, the complex nature can be eliminated by algebraic 
manipulation. In particular, the four nonzero eigenvalues 
are of the form 



( 151 ) 

The associated eigenvectors are of the form 

*i - ( ± g 2 ± i f 2 ) and <fi i = ( ± j 2 ± i * 2 ) • 

( 152 ) 

The nonsubscript i being I , not an index. Expanding, 

*2 * ^ g 2 + ^ ^ 2 ^ ' *1 = *2 = ^ g 2 “ * 

*4 = ( " g 2 " 1 f 2 ) ; *3 = = ( _g 2 + 1 f 2 ) 

* 2 " (j 2 + 1 k 2 } ? *1 = ~2 = (j 2 ’ 1 k 2 > 

0 4 = (-j 2 - i k 2 ) ; 0 3 - "0^ = (-j 2 + i k 2 ) 

( 153 ) 
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The constants assume the form 


w 2 = c 2 + 1 d 2 ; w i - C 1 + 1 d l 


w 4 = c 4 + i d 4 ; W 3 = c 3 + i d 3 

(154) 

Starting with the equation for the displacement w* (x) , eq. 
150, the above eigenvalue and eigenvector definitions are 
substituted. The displacement w* (x) then assumes the form 

(a - i b ) x 

w * (x) - ( c x + i d x ) e * 2 


( a + i b ) x 

+ ( c 2 + i <^ 2 ) e 


( -a + i b ) x 

+ < c 3 + 1 d 3 } e 


( -a - i b ) x 

+ ( c 4 + i d 4 ) e 


+ a U- + p V, + T 2 - . (155) 


Note that a and |3 are real . The real and imaginary parts 
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of the exponential terms are separated as follows: 


e (x + i y) * e x ( cos y + i sin y ) , (156) 


ia x 

e = cosh a,x ± sinh a_x , 


(157) 


and these are substituted into eq. 155. Since w* is real, 
the imaginary part of the right hand side of eq. 155 must 
equal zero, therefore, 


c 


2 


c 


1 ' 


c 


4 ' 




(158) 


Substitution and regrouping of this information into odd 
and even functions of x reveals 


w* (x) = 2 ( c 2 + c 4 ) cos b 2 x cosh a 2 x 


2 ( d 2 + d 4 ) sin b 2 x sinh a 2 x 


even 


+ a u 6 + * V s + N e 




« 
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odd —\ 


c. ) cos b n x sinh a n x 
4 2 2 


+ 2 ( C 2 - 

+ 2 ( d 4 - d 2 ) sin b 2 x cosh a 2 x . (159) 


Separation into odd and even functions of x was done 
because the axial compression, torsion, and temperature 
change problems all lead to responses that are 

characterized by w* (x) being an even function of x. 
Therefore, the odd part of w* (x) must be zero. For this to 
be true 


c„ = 



(160) 


Hence, the w* (x) displacement equation reduces to 


w* (x) - 2 C cos b 2 x cosh a 2 x 


2 D sin b 2 x sinh a 2 x + a U 


+ 9 v 6 


+ 



(161) 


where 



C - c 2 + C 4 
D = d 2 + d 4 . 


Expressions for 

u* (x) and 

V (x) 

can be 

derived 

similarly 

using the 

information 

found 

from the 

V (x) 

equation. 

particularly 

eqs. 158, 

160, 

and 162. 

These 


displacements are found to be 


u* (x) =* 2 C ( g 2 cos b 2 x sinh a 2 x - 


f 2 sin b 2 x cosh a 2 x ) 


2 D ( f 2 cos b 2 x sinh a 2 x + 


g 2 sin b 2 x cosh a 2 x ) + U g x ; (163) 


and 


v* (x) = 2 C ( j 2 cos b 2 x sinh a 2 x - 


k 2 sin b 2 x cosh a 2 x ) 


- 2 D ( k 2 cos b 2 x sinh a 2 x + 
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j 2 sin fc> 2 x cosh a 2 x ) + V g x . (164) 

By choosing w* (x) to be even in x, u # (x) and v* (x) are odd 
in x. This was expected. The above solutions for u e (x) , 
v* (x) , and w* (x) are valid for any axisymmetric response 
that is characterized by w* (x) being an even function of x. 
The boundary conditions for each of the three problems lead 
to specific solutions based on the above general forms. In 
the next three sections the boundary conditions for the 
three cases considered here are presented. 



Two sets of variationally consistent boundary 
conditions for axial compression of the cylinder will be 
considered. The first set of boundary conditions is 


• 

u* 

( L/2 ) 

= -0.5 t ; 

(165) 


V * 

( L/2 ) 

= 0 ; 

(166) 

• 






w* 

( L/2 ) 

- 0 

(167) 

• 

dw* 

( L/2 ) 





dx 

= 0 

(168) 
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The second set of boundary conditions is 


u* ( L/2 ) = -0.5 t ; 


( N xQ ( L/2 ) 

V ( L/2 ) = 

dW f L/2 ) 
dx 


+ 


M x0 1 L / 2 > 
R 


0 


0 . 


( 169 ) 


0 ; (170) 


(171) 


(172) 


It is important to note that only the conditions at one end 
of the cylinder are considered. By the symmetry conditions % 

on u* (x) , V (x) , and V (x) , conditions on the other end are 
automatically specified. 

The first set of conditions corresponds to a cylinder % 

that is restrained on each end against circumferential 
displacement, radial displacement, and rotation in the x 
direction. The ends of the cylinder are compressed toward 4 

each other with a displacement equal to a laminate 
thickness, t. The second set of conditions corresponds to 
a cylinder that is free to move circumferentially at each • 

end but is restrained against radial displacements and 
rotation in the x direction. The ends of the cylinder are 
again compressed toward each other. This compression ^ 
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condition is considered for the case where temperature 
effects are ignored and so N 0 T = 0 in eq. 117. Herein 
the above two sets of boundary conditions for compression 
will be referred to as boundary condition 1 and boundary 
condition 2, respectively. They will be abbreviated be 1 
and be 2. These two sets of boundary conditions are used 
in conjunction with eqs. 161-164 to solve for the four 
unknowns, C, D, U 6 , and V 6 . After values for the unknowns 
are found, values of the displacements and force and moment 
resultants along the cylinder can be calculated. 

3 . 8 The Boundary Conditions for Torsion 

Two sets of variationally consistent boundary 
conditions for twisting the cylinder are considered. The 
first set is 


U- ( L/2 ) 

= 0 

(173) 

V ( L/2 ) 

= 0.5 t ; 

(174) 

V ( L/2 ) 

- 0 

(175) 

dw° ( L/2 ) 

_ r\ 

(176) 

dx 

— u 
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The second set of boundary conditions is 


N x ( L/2 ) = 0 ; 

v* ( L/2 ) - 0.5 t ; 

V ( L/2 ) = 0 

dW ( L/2 ) = 0 

dx 


(177) 

(178) 

(179) 

(180) 


The first set of boundary conditions correspond to the case 
of a cylinder with its ends restrained against axial 
displacement, radial displacement, and rotation in the 
axial direction, and with a circumferential displacement of 
one end of the cylinder relative to the other equal to one 
laminate thickness. The second set is similar to the first 
except that the cylinder is free to expand or contract in 
the axial direction. These conditions are denoted as be 1 
and be 2. Here also the twisting is considered independent 
of temperature effects and so Nq^ = 0. Similar to 
compression, these two sets of boundary conditions are used 
to solve for the four unknowns C, D, Ug, Vg and then to 
calculate the values of the displacements and force and 
moment resultants along the cylinder. 
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3 . 9 The Boundary Conditions for Thermal Expansion 


There is only one set of boundary conditions that 
correspond to the case of free thermal expansion of the 
cylinder. They correspond to the ends of the cylinder 
being completely free of any applied loads. These are 
given by the variationally consistent boundary conditions 

N x ( L/2 ) - 0 ; (181) 

M x ( L/2 ) - 0 ; (182) 


dM x ( L/2 ) 
dx 


0 


(183) 


N xe ( L/2 ) - 0 . (184) 

Substituting the expressions for the displacements egs. 
161-164 into the force and moment resultant expressions, 
eqs. 109, 111, and 112; and substituting these results into 
the above equations leads to four equations from which to 
solve for the four unknowns, C, D, U 6 , and V 6 . 

The next chapter outlines the three FORTRAN programs 
used to solve for the four unknowns, C, D, U 6 , and V 6 . 



Chapter 4 - Creation of the Software 


The preceding analyses were implemented into three 
FORTRAN computer programs, one for each of the three types 
of loading conditions. All three programs have the same 
logic, only the boundary conditions and some input 
information are different. Figure 6 illustrates a 
simplified flowchart of the programming logic. The 
programs were written for minimal user input but with as 
much generality as possible. Therefore, only the cylinder 
geometry, which includes the length and mean radius, and 
the material properties, which include Young's moduli in 
the x and y directions, shear modulus, Poisson's ratio, the 
thickness of each layer of the material, the number of 
layers, and the orientation of the fibers of each layer, 
are input. For thermal expansion, the two coefficients of 
thermal expansion, a x and ay, and the change in temperature 
are also input. After input, the reduced and transformed 
reduced stiffness matrices are calculated. The calculation 
of the transformed reduced stiffness matrix enables the ABD 
matrix to be calculated. The input and calculated values 
are then written to a data file. Numerical values for C 8 , 
Cg, and C 4 of eqs. 13 0-132 are calculated using the above 
information. With numerical values for these variables, 
the roots or eigenvalues are found by solving eq. 133. 
From the eigenvalues, the eigenvectors are found directly 
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using Cramer's rule and a closed-form expression for the 
determinant of a 4 x 4 matrix. With the eigenvalues and 
eigenvectors known, the unknown constants can be found for 
the particular set of boundary conditions. The four 
unknowns are also found directly using Cramer's rule and a 
closed- form expression for the determinant of a matrix. 
Using the three equations for displacement, eqs. 161-164, 
and the six stress resultant equations, eqs. 109-113, 100 
data points are calculated along the cylinder from the 
center at x * 0 to the end at x =» + L/2. Each set of data 
points is written to its own data file for each of the 
displacements and stress resultants of each boundary 
condition. These files are used for plotting the results. 

The listing and/or a copy on a floppy disk for each 
of the three programs is available upon request. 



Chapter 5 - Results 


To determine the character of the responses that are a 
result of the elastic couplings and directional dependence 
of fiber-reinforced materials, several specific examples of 
cylinders are studied for the three loading cases. The 
first example discussed could be considered somewhat 
unusual. However, it was chosen to determine the types of 
responses that can be exhibited by composite cylinders. 
The cylinder studied is a two layer cylinder with a 
stacking arrangement of (23/-60) T . In the stacking 
sequence notation the number denotes the angle, in degrees, 
of the fibers relative to the axial direction of the 
cylinder. The first entry in the notation denotes the 
orientation of the inner layer, and subsequent entries 
correspond to layers further away from the inner radius. 
Clearly, such a cylinder is highly coupled in that it is 
both unbalanced and unsymmetric, i. e., A 16 f 0, A 2 g ? 0, 
and there are entries in the B^j matrix. The other 
examples are more realistic and do not deviate too far from 
current design philosophies. However, these latter 
examples exhibit some characteristics not usually found in 
current designs. 

The material properties of each layer are assumed to 
represent AS4/3501 graphite-epoxy. This material cures at 
350 * F and has the following elastic properties: 
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20. E6 psi 


/ 


E 


x 


Ey = 1.3 E6 psi 


i 


E = 1.03 E6 psi 

s 

= ' 3 


In each case, the cylinder geometry is: 

mean radius R = 10.0 inches , 

length of cylinder L = 30.0 inches 


To illustrate the response of the various cylinders, the 
displacement and stress resultant responses of each 
cylinder as a function of axial distance along the cylinder 
is illustrated in a series of figures. The displacement 
response of the cylinders is normalized by the laminate 
thickness and axial distance is normalized by the length of 
the cylinders. Since the responses of the right and left 
halves of the cylinder are either even or odd functions of 
axial distance, the responses for just the right half of 
the cylinder are shown. 
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5 . 1 Response of the (23/-601 T Cylinder 


5.1.1 Compression 

The values of the various entries in the ABO matrix 
for the ( 23/-60 ) ij laminate are: 


A matrix 


A 11 

- 

.90107 E5 lb/ in 

A 12 

- 

.30097 E5 lb/in 

A 16 

* 

.15812 E5 lb/ in 

A 22 

- 

.71799 E5 lb/ in 

A 26 

* 

-.22710 E5 lb/in 

A 66 


.36474 E5 lb/ in 

B matrix 


B 11 


-152.45 lb , 

B 12 

s 

12.006 lb , 

B 16 

as 

-96.633 lb , 

B 22 

- 

128.44 lb , 

B 26 

- 

-89.747 lb , 


(185) 
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12.006 lb 


(186) 


and the D matrix 


D 11 

= 

.75089 

lb-in 

°12 

25 

.25081 

lb-in 

°16 

= 

.13177 

lb-in 

°22 

* 

.59832 

lb-in 

°26 

- 

-.18925 

lb-in 

D 66 

- 

.30395 

lb-in 


9 


9 


9 


9 


(187) 


As can be seen, there are no zero entries in the matrix, 
this laminate construction exhibiting full coupling. 

The axial displacement response of the (23/-60) T 
cylinder for the case of axial compression with the two 
boundary conditions is illustrated in fig. 7. The response 
is linear with x, a characteristic which is due to the 
dominance of the Ug term of eq. 163, and starts at zero at 
mid-cylinder and goes to minus one-half at the right end. 
Recall, the compression condition is achieved by displacing 
one end of the cylinder relative to the other by one 
laminate thickness. Therefore, the right end moves 
one-half a laminate thickness in the negative x direction. 
As can be seen, the axial displacement response of the 


74 


UO DISPLACEMENT 


(23/— 60)T, UO VS X 



Figure 7 (23/-60) T , u* (x) vs x, Compression 


75 



cylinder is independent of boundary condition. On the 
other hand, the axial force resultant, illustrated in fig. 
8, is strongly dependent on the boundary condition. The 
difference in the magnitudes of the axial force resultant 
in fig. 8 is due to differences in the circumferential 
boundary conditions. For metallic and balanced symmetric 
composite cylinders such a dependence simply does not 
exist. Circumferential displacements do not strongly 
influence, if they influence at all, the axial force 
characteristics. The highly coupled nature of the 

responses is a characteristic that is unique to 

fiber-reinforced material. This can be clearly seen by 
examining eq. 109, the expression for the axial force. The 
third term in the equation involves gradients of the 
circumferential displacement. In metallic and balanced 
symmetric composite cylinders, A 16 , B 16 , and B 1;1 are all 

zero and so the axial force involves only the axial 
displacement gradient, through the first term in that 
equation, and the radial displacement, through the second 
term. The second term, as mentioned previously, represents 
the standard coupling of radial and axial responses due to 
Poisson effects. The results of fig. 8 indicate that the 
cylinder restrained against twist, be 1, requires a larger 
axial compressive force to move the ends of the cylinder 
toward each other than the cylinder that is free to twist. 
The difference in the two force levels is a direct result 



(Ib/ln) 





(23/— 60)T, NX LOAD VS X 


COMPRESSION 




Figure 8 (23/-60) T , N x vs x, Compression 
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of this twist term, specifically the A 16 term. 

To further illustrate the coupling of axial and 
circumferential responses, the circumferential displacement 
response of the cylinder for the two boundary conditions is 
shown in fig. 9. As expected, the circumferential 
displacement is strongly influenced by the type of boundary 
condition. Boundary condition 2 imposes no constraint on 
circumferential displacement at the ends of the cylinder 
and so the cylinder twists as it is compressed axially. As 
was seen in fig. 8, this freedom to twist reduces the axial 
force required to compress the cylinder. The twist 
response is a characteristic of V 6 in eq. 164. Boundary 
condition 1, on the other hand, restrains twist motion. 
The restraint is not complete because there is actually 
twisting of the cylinder between the ends and center of the 
cylinder. The twist displacements are not discernible on 
the scale of fig. 9. To complement fig. 9, the shear force 
resultant N x0 for the two boundary conditions is 
illustrated in fig. 10. This figure, like the previous 
three, typifies the coupling between axial and 
circumferential responses. To compress the 
circumferentially restrained cylinder, be 1, requires a 
shear force resultant. The expression for N x0 , eq. Ill, 
clearly illustrates this. The constant A 16 multiplies the 
axial strain, the axial strain being the dominant response 
for this case of compression. 


For the case 



VO DISPLACEMENT 


• (23/— 60)T, VO VS X 



Figure 9 (23/-60)<j», v* (x) vs x, Compression 
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Figure 10 (23/-60) T , N x0 vs x, Compression 9 
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of a restrained cylinder the third term, the term normally 
associated with twist response, is zero, and so the 
majority of the contribution to N x @ is due to the axial 
strain term. For the unrestrained cylinder, the twist 
response is not zero and the third term cancels the first 
term in N xe . Note that A^g was responsible for the 
coupling in the expression for N x , eq. 109. 

The radial expansion of the cylinder, i.e., the w° 
displacement response, for the case of axial compression is 
shown in fig. 11. Due to Poisson effects, the axial 
compression of any type of cylinder, metallic or composite, 
results in radial expansion. Here, however, the various 
elastic couplings cause the cylinder with no twisting 
restraint to expand more. This results in a somewhat 
larger displacement gradient in the boundary layer at the 
supported end of the cylinder. The displacement gradient 
in the boundary layer influences the circumferential force 
resultant N e , as shown in fig. 12. Over the majority of 
the cylinder there is no change in the value of the radial 
displacement. The value of the circumferential force 
resultant is zero there. However, at the end of the 
cylinder, the radial displacement changes rapidly causing a 
circumferential force resultant there. For the cylinder 
with no twist restraint, this induced stress resultant is 


larger in 

magnitude 

than 

when 

the 

twist restraint 

is 

present. 

Examining 

eq. 

110 

it 

is seen that 

the 






Figure 12 


(23/-60)«p, Nq vs x. Compression 
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circumferential force resultant and the circumferential 
displacement are coupled through A 26 and B 26 in the third 
term. These coefficients are zero in metallic and balanced 
symmetric composite cylinders. For the circumferentially 
restrained cylinder dv'/dx is zero and the third term has 
no contribution. In the unrestrained cylinder, there is a 
gradient in v* and this contributes to the circumferential 
force resultant and, as noted above, causes it to be larger 
than when twist is restrained. 

The bending moment resultants M x and M e are 
illustrated in figs. 13 and 14, respectively. Referring to 
egs. 112 and 113, for symmetrically laminated composite 
structures all coefficients in these equations except the 
D^j are zero. For some symmetrically laminated composite 
structures D^g and D 2 g are small, and for metallic 
structures they are zero. Therefore, the bending moments 
depend strongly on the radial displacement only through the 
axial curvature of the radial displacement, the fourth 
term. The circumferential and axial displacements have no 
influence on the moments. Clearly this is not the case for 
general composite laminates. The axial and radial 
displacements enter through the B^j terms and influence the 
moments. In addition, the moments M x and M e are strongly 
influenced by the presence or lack of circumferential 
displacement, i.e., dv*/dx. For both M x and M e , the first 
three terms in the moment expression are relatively 


MX MOMENT (ln-lb/tn) 


0.16 


(23/— 60)T, MX MOMENT VS X 


COMPRESSION 



X 


Figure 13 


(23/-60) M x vs x. Compression 
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(23/— 60)T, 



MOMENT VS X 



Me vs x, Compression 



independent of x. Over the majority of the cylinder the 
curvature is zero and so the fourth terms are zero. 
However, near the supported ends, the curvature becomes 
large and the fourth terms dominate both moment 
expressions. This dominance is confined to a small region 
near the end and thus constitutes a boundary layer there. 
For the case of Me and be 2 , there is no moment resultant 
except for in the boundary layer. 


5.1.2 Torsion 

As with compression, the twisting of one end of the 
(23/-60) <p cylinder relative to the other results in several 
unexpected responses. The circumferential displacement as 
a function of axial distance from the center of the 
cylinder is illustrated in fig. 15. This displacement is 
linear with distance, and is the counterpart to fig. 7 for 
the axial compression case. The results for both the 
axially restrained case, be 1, and the axially free case, 
be 2, are identical, the response being determined by v 6 in 
eq. 164. The shear force resultant, N xe , required to 
produce the twist is shown in fig. 16. There is 
considerable difference in the shear force required for the 
two boundary conditions. The presence of an axial 
restraint considerably stiffens the cylinder in torsion, a 



(23/— 60)T, VO VS X 


TORSION 




NXTHETA LOAD (Ib/ln) 


12.00 


(23/— 60)T, NXTHETA LOAD VS X 

TORSION 



X 


Figure 16 


(23/-60) <j, N xe vs x > Torsion 
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direct result of the elastic coupling within the cylinder. 
Examination of the expression for the shear force 
resultant, eq. Ill, indicates that the axial displacement, 
or lack of it, enters the problem through the coefficients 
A]_g and B^g. The character of the axial displacement for 
the two conditions is shown in fig. 17, the axially 
unrestrained case exhibiting a linear variation of axial 
displacement with length and the axially restrained 
cylinder showing very little variation with length. As 
with compression, this linearity in the unrestrained case 
is due to the dominance of the Ug term of eq. 163. The 
presence or lack of axial displacement causes, through A 16 
and Bj^g, the difference in shear force resultant values. 
The axial-twist coupling is further illustrated by 
examining the difference in the axial force resultant 
caused by the twist for the two cases. It is indicated in 
fig. 18 that to restrain the axial motion of the cylinder 
during twist, be 1, an axial force the same order of 
magnitude as the force required to compress the torsionally 
restrained cylinder, fig. 8, is necessary. It should be 
noted that neither N x nor N x0 exhibit boundary layers in 
either of the two loading cases studied so far. This is 
because the governing differential equations, eqs. 95 and 
96, indicate the two resultants cannot vary with x. They 
must assume the boundary value along the complete length of 
the cylinder. 


i 




90 


I 


UO DISPLACEMENT 


(23/— 60)T, UO VS X 


TORSION 



Figure 17 (23/-60) T , u* (x) vs x, Torsion 
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(23/— 60)T, NX LOAD VS X 
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Figure 18 (23/-60) T , N x vs x, Torsion 



For the case of twist, perhaps the most unusual 
effect is the coupling between twist and radial 
displacement, fig. 19. Because the laminate is 
unsymmetric, these two displacements are coupled, 
particularly through the B 16 term in eqs. Ill and 112. It 
is shown in the figure that twisting of the cylinder causes 
its radius to increase. Due the boundary condition 
restraining radial displacement at the ends of the 
cylinder, there is a boundary region in the radial 
displacement at the ends. Figure 19 looks very similar to 
fig. 11, the radial displacement caused by the axial 
compression. The character of fig. 11 would be much the 
same for the case of a metal or balanced symmetric 
composite cylinder. However, there would be no radial 
displacement due to twist for a metal cylinder and so the 
results of fig. 19 are unique to composites, and in 
particular, unsymmetric composites. Furthermore, there is 
a sign sensitivity to the coupling in that if the cylinder 
was twisted the other direction, the radial displacements 
would be negative. In addition, the sign of B^g can be 
reversed by reversing fiber angles, i.e., fabricating a 
(-23/60) <p laminate. Then the radial displacement would 
also decrease. This radial displacement-twist coupling is 
felt to be the most significant effect found to this point 
in the study. 

The character of the moment resultant normally 



(23/— 60)T, WO VS X 


TORSION 



Figure 19 (23/-60) T , V (x) vs x, Torsion 
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associated with the radial displacement, namely M x , is 
illustrated in fig. 20. Due to bending in the axial 
direction, there is an increased moment at the end of the 
cylinder. It is interesting to note that even though there 
is no curvature in the axial direction over the majority of 
cylinder length, i.e., w* is a constant over much of the 
cylinder, there is a moment generated. This is because 
with composites, moments are associated with displacement 
gradients other than just curvatures, in particular, 
gradients of inplane displacements, and with displacements 
themselves. This has been discussed before and can be seen 
by examining egs. 112 and 113. With metals and balanced 
symmetric laminates, it is only curvatures that are 
associated with moments. 

The character of the other stress resultants 
associated with twisting the cylinder are illustrated in 
figs. 21 and 22. As with the axial compression case, the 
circumferential force resultant Nq increases at the end of 
the cylinder, and the moment resultant M@ is small and 
constant over the majority of the cylinder, exhibiting an 
increased value in the boundary layer region. 


5.1.3 Temperature Change 


Due to the elevated temperature processing associated 
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with most polymer-matrix composites, the thermally-induced 
deformations and stress resultants, from which 
thermally-induced stresses can be determined, should be 
computed to determine their influence on the initial shape 
of the cylinder, or the initial stress state. For a 
linear theory, which is used here, these deformations and 
stress resultants superpose on the response due to axial 
compression or twist. If deemed important, the 
thermally- induced responses can be computed in a separate 
step and the results added to the results for the 
mechanically- induced response. For a nonlinear theory, the 
thermally- induced deformations and stress resultants are an 
important first step and should be the starting point in 
any analysis directed at understanding the response due to 
loads. If the loads are being applied at a temperature 
lower than the processing temperature, deformations and 
stresses are present before the load is applied. For a 
nonlinear theory, these thermally-induced responses do not 
simply add to the nonlinear response. They are a more 
integral part of the solution. Including or not including 
the thermally-induced responses as the initial starting 
condition from which to apply the mechanical load may lead 
to entirely different predictions for the response due to 
the load. In this study the intent is to simply examine 
the thermally- induced response itself, without adding it to 
the mechanically- induced results, or implying what the 



thermal response may mean as regards nonlinear effects of 
mechanical loads. 

The processing temperature of AS4/3501 is 350 *F and 
if it is assumed that the cylinder is at a service 
temperature of 70 *F, this represents a -280 *F temperature 
change. The axial displacement response of the cylinder as 
a result of this temperature change is shown in fig. 23. 
The axial displacement of the right-half of the cylinder is 
negative everywhere, meaning the cylinder length contracts 
due to the temperature change. The displacement response 
is linear with position except at the extreme end, where 
there is a slight perturbation from the linear response. 
The circumferential displacement as a function of axial 
position is shown in fig. 24. This response is also linear 
with axial position except for a slight perturbation at the 
end. The reason for the perturbation is quite clear when 
the radial displacement response, fig. 25, is examined. 
Away from the ends, the cylinder contracts uniformly in the 
radial direction as the temperature is lowered from the 
processing temperature. However, at the end of the 
cylinder the radial displacement changes rapidly from this 
uniform state. As the end of the cylinder is approached 
the radius decreases more, and then suddenly increases. 
This change from the uniform state occurs over a very short 
axial distance, resulting in large displacement gradients 
and thus a boundary layer in this region. 
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Figure 25 


(23/-60)ij, w* (x) vs x, Thermal Expansion 
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There are thermally- induced loads associated with 
this deformed shape. Due to the traction- free boundary 
conditions, the axial force resultant, the shear force 
resultant, and the axial moment resultant, N x , N xe , and M x , 
respectively, must be zero at the end of the cylinder. The 
other resultants are not specified anywhere along the 
cylinder. As can be seen in fig. 26, the axial moment 
resultant along the length of the cylinder is uniform and 
nonzero over much of the length, increases slightly, then 
rapidly decreases to zero as the end of the cylinder is 
approached. If the (23/-60) T laminate was fabricated and 
processed as a flat plate and then cooled to 70 *F, the 
unsymmetric nature of the laminate would cause it to warp 
significantly and develop large curvatures [1,2,3]. 
However, the cylinder cannot develop curvatures in the 
circumferential direction nor in the axial direction except 
near the end where it is traction free. With the laminate 
naturally wanting to develop curvature but not being able 
to because of geometric constraints, force and moment 
resultants develop in the cylinder. However, with the end 
of the cylinder being traction free, three of the force and 
moment resultants must be zero there. If they are not zero 
away from the ends, they must ’adjust themselves' and 
vanish at the end. This is the case with M x . As a result 
an axial curvature develops in the radial displacement 
near the end of the cylinder. 
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Figure 26 (23/-60) T , M x vs x, Thermal Expansion 
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The circumferential force resultant N e , fig. 27, 
behaves in the opposite manner. This force resultant is 
zero over most of the length of the cylinder. Toward the 
end the resultant becomes negative and then rapidly becomes 
positive. Since there are none, this behavior is not 
associated with any imposed kinematic constraint at the end 
of the cylinder, as it is with the axial compression and 
twist cases. Rather, this behavior is due to the coupling 
of the stress resultants through the various displacement 
components, most notably the radial displacement. The 
circumferential moment resultant, Mq, is uniform over most 
of the cylinder and then varies rapidly near the end, as 
illustrated in fig. 28. The axial and shear force 
resultants are exactly zero over the entire length of the 
cylinder. As stated before, the governing equations 
stipulate that neither of these resultants vary with x, and 
since they are zero at the ends, they are zero everywhere 
along the length. 

In summary, it is clear the (23/-60) T cylinder 
exhibits a high degree of elastic coupling. Though it is a 
particular case, this cylinder typifies the fact that with 
composite materials, twist, extension, and bending can be 
coupled. The most significant finding is the coupling of 
radial displacement with twist displacement for the 
cylinder subjected to torsion. 
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Figure 27 


(23/-60) <p, N e vs x, Thermal Expansion 
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(23/— 60)T, MTHETA MOMENT VS X 


THERMAL EXPANSION 



Figure 28 (23/-60)^, M e vs x, Thermal Expansion 
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5 . 2 Response of a Ouasi-isotropic Cylinder and Two 

Variants 

The results for the (23/-60) T cylinder indicate the 
wide range of responses that are possible with general 
composite laminates. The coupling of the radial 
displacements with the axial and circumferential 

displacements, and the accompanying coupling of the force 
and moment resultants, are characteristics that can be 
tailored for particular applications. The (23/-60) T 
laminate is quite unorthodox when compared to current 
designs. Unfortunately, unorthodox designs can often be 
viewed with skepticism due simply to lack of familiarity. 
Positive features are often overlooked. Therefore, it is 
of value to examine slight variations in existing designs, 
and explore the characteristics of these designs. 

Specifically, a popular conventional design is the 
quasi-isotropic laminate. This laminate consists of equal 
percentages of +45*, -45*, 0*, and 90° laminae, the 
percentage being 25% for all four angles. In the plane of 
the laminate the effective Young's modulus, Poisson's 
ratio, and shear modulus are independent of direction, 
i.e., isotropic. However, the shear modulus is not related 
to Young's modulus and Poisson's ratio by the well known 
relation 




G 


E 

2 ( 1 + H ) 


(188) 


Because of the lack of a relation between these three 
properties, the term quasi-isotropic is used. A typical 16 
layer quasi-isotropic laminate is given by the stacking 

arrangement 

( ±45/0/90 ) 2s 

where the subscript s implies the laminate stacking 
arrangement is symmetric with respect to the laminate 

midplane. Thus both the lamina at the inner radius and the 
lamina at the outer radius are oriented at +45*. The 
laminae on either side of the mean radius are oriented at 
90*. Quasi-isotropic designs are popular because many 
metal-based design concepts are usable and because of the 
lower density of graphite/polymer composites, a simple 
substitution of a quasi-isotopic laminate for aluminum 

saves 25% weight. Accordingly, then, except for the 

presence of the D^g and D 2 g terms, a quasi-isotropic 
cylinder behaves very much like a metallic cylinder. 

Two variants of the above quasi-isotropic design are 

| 

i 

design 1 ( +45/-45/0 4 /+45/-45 2 /+45/90 4 /-45/+45 ) T 

1 .... (0° layers closer to inner radius) 
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design 2 ( +45/-45/90 4 /+45/-45 2 /+45/0 4 /-45/+45 ) T 

.... (90* layers closer to inner radius) 

These designs have the same number of layers as the above 
quasi-isotropic design and these layers have the same 
orientations as the quasi-isotropic design. In fact, all 
three laminates have the same A matrix, the matrix 
associated with inplane stiffness, and so all three 
laminates are quasi-isotopic in their inplane properties. 
However, both designs 1 and 2 are unsymmetric. Design 1 is 
distinguished by having the 0* layers grouped together and 
placed close to the inner radius. Design 2 has the 90* 
layers grouped and placed close to the inner radius. The 
placement of the 0* and 90* layer groups provides some 

eccentricity in the material properties. The eccentricity 
of design 1 is exactly opposite to the eccentricity of 
design 2. The motivation for considering eccentricity is 
as follows: It is well known that the axial buckling 

response of a cylinder is influenced by having axial 
stiffeners on the inside or outside of the cylinder. Among 
other things, the stiffeners provide eccentricity and can 
delay the onset of buckling. For a certain range of 
geometric parameters, stiffeners on the inside are more 

beneficial than stiffeners on the outside. For other 
ranges of parameters, stiffeners on the outside are most 

beneficial. Whether the stiffeners are on the inside or 


outside determines the sense of the eccentricity. The idea 
here is to provide the eccentricity by material coupling 
rather than by the addition of stiffeners. This may save 
weight and certainly simplifies fabrication. A complete 
understanding of whether this material coupling is 
beneficial cannot be determined until a buckling analysis 
is conducted. That is beyond the scope of this study but 
it is felt examining the linear response of unsymmetric 
cylinders is a worthwhile starting point. It is also 
worthwhile determining the differences between similar 
symmetric and unsymmetric designs. That is done in this 
section. 

The A matrix of all three designs is: 


A 11 = A 22 = • 69181 E6 Win 

A 12 - .19659 E6 lb/in , 

A 16 " A 26 ** ib/in i 

A 66 = .24761 E6 lb/in , (189) 


The B matrix of a symmetric laminate is of course zero and 
the B matrix of design 1 is: 


B 11 “ B 22 


-7524. lb 


B 12 “ B 16 B 26 B 66 


0. lb 


§ 


(!90) 


while the B matrix of design 2 is: 


B 11 “ - B 22 = 7524 ' lb ' 

B 12 * B 16 * B 26 - B 66 - °- lb - U 9 U 

The D matrices of the symmetric design, and design 1 and 
design 2 are, respectively, 


'll 

_ 

364.14 

lb-in , 

°12 = 

137.89 

lb-in , 

16 

=* 

23.513 

lb-in , 

°22 = 

307.71 

lb-in , 

*26 

S 

23.513 

lb-in , 

°66 

165.10 

lb-in , 







(192) 

'll 

= 

352.44 

lb-in , 

°12 = 

121.37 

lb-in , 

16 

s 

18.810 

lb-in , 

°22 = 

352.44 

lb-in , 

26 

= 

18.810 

lb-in , 

°66 

148.58 

lb-in 


(193) 


As can be seen, the major difference between designs 1 and 
2 is in the B matrix. .. .the signs of B^ being reversed. 
The sign reversal on B 1;l is due to a reversal of the 
material eccentricity. It should be noted that the B 
matrix for either case is sparsely populated, and B 2 2 = 
- B^. At this point it is interesting to compare the ABD 


matrices for these three laminates, eqs. 189-193 with the 
ABD matrix of the (23/-60) T laminate, eqs. 185-187. The 
two main differences are: (1) For the (23/-60) T laminate, 
the B matrix is completely full and contains no zero 
elements. For designs 1 and 2, all elements in the B 
matrix except 2 are zero, and these two are of equal 
magnitude. (2) For the symmetric design and design 1 and 
2, Aig and A 2 6 are zero, for the (23/-60) T case they are 
not. The lack of the A^ 6 and A 26 terms means extension and 
twist cannot directly couple. The D matrices are full for 
all cases (a characteristic of most laminates) . 

The next sections describe the response of the three 
designs to the compression, torsion, and a temperature 
change. Any differences in the response among the three 
designs are noted. 


5. 3 Response to Compression 

The axial displacement for the three cylinders as a 
function of length along the cylinder is illustrated in 
fig. 29. The boundary conditions considered are boundary 
condition 1, eqs. 165-168, namely the cylinder restrained 
against circumferential displacement at the ends. As 
expected, the response is quite linear with x with the U 6 
term being dominant, and the three cylinders exhibit nearly 






I 


identical response. Though not visible on the scale of the 
figure, the differences that do occur are at the end. The 
axial force resultant required to compress the cylinder one 
laminate thickness is shown in fig. 30. The three cases 
are shown and these forces are identical. The forces are 
constant with x, because as mentioned earlier, the 
differential equations require this. 

Slight differences in the response are apparent when 
the radial displacements are studied. This is illustrated 
fig. 31. All three cylinders have identical radial 
displacement away from the end and exhibit a boundary layer 
at the end. The unsymmetric design 1, the design with the 
0* laminae closest to the inner radius, results in the 
least amount of radial displacement variation within the 
boundary layer. Design 2 exhibits the most, and the 
symmetric design is between the two. Overall, design 1 is 
more like the symmetric design than design 2 is. If it is 
to have an effect at all, one might expect the smaller 
gradient and smaller reverse curvature of design 1 within 
the boundary layer to delay the onset of buckling, whereas 
the larger gradient and larger reverse curvature of design 
2 might aggravate buckling. 

The axial moment resultants that accompany the three 
radial displacement characteristics are shown in fig. 32. 
Away from the ends there is no axial moment generated in 
the symmetric cylinder. The only moment generated is in 
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Figure 30 Three Designs, N x vs x. Compression 
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Figure 31 Three Designs, w* (x) vs x, Compression 
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THREE DESIGNS, MX MOMENT VS X 



Figure 32 Three Designs, M x vs x, Compression 





the boundary layer. However, the axial compression of both 
unsymmetric designs generates moments away from the ends, 
the moments being exactly opposite. The moment away from 
the end is altered as the end of the cylinder is 
approached. Except for slight differences at the ends the 
three moment characteristics are essentially a positive or 
negative vertical shift of symmetric case. 

The circumferential force resultants for the three 
cases, fig. 33, differ very little. On the otherhand, the 
circumferential moment resultants shown in fig. 34 are 
quite different. In addition to having a nonzero value 
over much of the cylinder length, unsymmetric design 1 
shows a large change in the circumferential moment near the 
end. Design 2 also exhibits rapid changes in the value of 
the circumferential moment near the end. It is interesting 
that all three designs converge to the same value of 
circumferential moment at the end of the cylinder, 
independent of the value of the moment at midcylinder, and 
independent of the local character of the variation within 
the boundary layer. 

Finally, the shear force resultant that accompanies 
axial compression is not zero but is very small ranging 
from 10" 10 to 10“ 9 . The presence of the term in eq. 
Ill is responsible for this. 
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Figure 34 Hg vs x, Compression 
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5 . 4 Response to Twisting 


The twisting responses of the three cylinders to 
boundary condition 1, eqs. 173-176, cylinders restrained 
against axial motion, are illustrated in figs. 35-38. The 
twisting of the three different cylinders does not produce 
as much response as the twisting of the (23/-60 )t cylinder. 
Most notably absent for all three designs is any radial 
displacement. Reviewing eg. 161, for these cases the 

unknowns C and D are exactly zero as is the constant {5 . 
Equation 145 indicates that 0 * 0 because A 2 $ = B 26 = 0. 
This makes the radial displacement exactly zero for all 

three designs. As seen in fig. 35, the circumferential 
displacement of all three cylinders is very similar, the 
variation of displacement being quite linear with position 
along the cylinder. As stated above the unknowns C and D 
are exactly zero. Since these unknowns are zero, the only 
term in eq. 164 to produce a response is the Vg term. 
Therefore, the response is linear from this term. The 
shear force resultants, N x q, fig. 36, are also quite 
similar for the three cylinders. The similar results are 
to be expected if the expression for N xe , eq. Ill, is 

examined. The material constants A 16 , A 26 , B 16 , and B 26 


are zero 

for 

all three 

cylinders 

and 

the shear force 

depends 

only 

on a 66' 

D 66 , and 

the 

circumferential 


displacement gradient. Since A 66 , D 66 , and the gradient 
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Figure 36 Three Designs, N x0 vs x, Torsion 
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are similar for the three cylinders, the shear force 
resultants are similar. As might be expected, because of 
its metal-like behavior, there is no axial or radial 
displacement response for the quasi-isotropic design. 
Recalling eqs. 161 and 163, these displacements are exactly 
zero which is due to the unknowns C, D, Ug being exactly 
zero for all three designs. As a result there is no axial 
or circumferential force resultant. The axial and 
circumferential bending moments, figs. 37 and 38, although 
small are not exactly zero. The reason can be seen by 
examining eqs. 112 and 113. The only nonzero terms are the 
D 16 and D 26 terms multiplied by the circumferential 
gradient. As mentioned previously and as seen in fig. 35, 
this gradient for the three designs is very similar and is 
the primary response that causes the bending moment. The 
difference between the symmetric and unsymmetric designs 1 
and 2 is in the D^g and D 2 g terms. This material property 
term differs slightly between the symmetric and unsymmetric 
cases, causing slight differences in the bending moments. 


5 . 5 Response to a Temperature Change 

Again, since graphite-epoxy is processed at an 
elevated temperature, the effects of this elevated 
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temperature on the displacements and stress resultants when 
the cylinder is at a lower temperature are of interest. 
Assuming the same -280 ’F temperature change as with the 
(23/-60) ip laminate, the response of the three cylinders is 
shown in figs. 39-43. The axial contraction due to the 
negative temperature change is illustrated in fig. 39. The 
axial contraction response of the three cylinders is 
practically identical, except at the ends where there are 
slight differences. These differences are much more 
significant when the radial displacement is examined. This 
is shown in fig. 40. Overall, the three cylinders 
contract radially. With the symmetric design, the 
contraction is uniform along the entire length of the 
cylinder. The two unsymmetric designs exhibit behavior 
different than the symmetric design, and different from 
each other. The radial displacement of design 1 actually 
becomes positive right at the end of the cylinder. The 
transition from uniform contraction to expansion takes 
place in a short boundary layer region. In fact, within 
the boundary layer, the radial contraction becomes more 
pronounced before suddenly changing to a radial expansion, 
with design 2 on the otherhand, the uniform contraction 
decreases as the end of the cylinder is approached, before 
suddenly becoming greater. The behavior of the unsymmetric 
laminates is due solely to the presence of the and B 2 2 


terms. 
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Figure 40 Three Designs, V (x) vs x, Thermal Expansion 
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The axial moments that accompanies the radial 
deflections are illustrated in fig. 41. No moment is 
induced in the symmetric laminate, whereas both unsymmetric 
laminates experience an axial moment. The moments for both 
cases are uniform over most of the cylinder length and 
suddenly go to zero in the boundary layer. With the 
symmetric cylinder there are no thermally- induced 
circumferential displacements, nor is there a shear force 
resultant. For the unsymmetric cylinders the 
circumferential responses are not exactly zero but are on 
the order of 10" 8 and 10“ 9 , respectively. The resultants 
in the circumferential direction consist of both forces and 
moments for the two unsymmetric cases, figs. 42 and 43. 
For the symmetric case, there are no circumferential stress 
resultants. For the two unsymmetric cases the 
circumferential force resultants are zero over most of the 
cylinder, and then rapidly become nonzero, changing signs 
once near the end of the cylinder. The moment resultants 
are of opposite magnitude over most of the cylinder, then 
change rapidly to dissimilar values at the ends. 

In summary, unsymmetric designs 1 and 2, though 
showing more coupling than the symmetric case, are not as 
highly coupled as the (23/-60)^> cylinder. This indicates 
that the degree of coupling can be controlled, and because 
of the similarity of some of the responses of the 
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Figure 41 Three Designs, M x vs x, Thermal Expansion 


133 



Figure 42 
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symmetric case and the two unsymmetric cases, designs 1 and 
2, certain responses obtainable with symmetric laminates 
can be retained. The most significant finding for the 
three cases just discussed is the difference in the radial 
displacements in the boundary layer when the cylinders are 
subjected to axial compression, i.e., fig. 31. 


Chapter 6 Conclusions and Recommended Future Research 


Presented has been the development of equations 
governing the axisymmetric response of fiber-reinforced 
composite cylinders to compressive, torsional, and thermal 
loads. The composite cylinders have been assumed to be 
general, exhibiting the highest degree of elastic couplings 
possible. The governing equations were developed by 
assuming Sanders' kinematic hypothesis was valid, and by 
using variational methods. Several sets of boundary 
conditions were considered, depending on the problem. 
Three FORTRAN programs were developed to implement a 
closed-form solution to the governing equations. Numerical 
results obtained from the programs for several example 
problems were presented and discussed. Based on examining 
the governing equations and the numerical results, several 
important conclusions can be drawn. These are as follows: 
For the most highly coupled case (the (23/-60)rp 
laminate) , a response totally unlike the response of metal 
cylinders or more conventional composite cylinders is 
exhibited. Compressive loading of a cylinder causes 
torsional stress resultants or circumferential 
displacements depending on the boundary conditions, and 
torsion loading causes axial stress resultants or 
displacements, again depending on boundary conditions. 
Furthermore, compressive loading of a cylinder causes 
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radial displacements beyond the normal displacements 
associated with Poisson effects. Most significant, 
torsional loading of a cylinder causes radial 
displacements. . . .an unusual effect. Boundary layer 
response is exhibited in most cases, the degree of response 
in the boundary layer depending on the particular case. 
Typically the boundary layer is characterized by large 
curvatures and sign reversals of the responses. It does 
not appear that the boundary layer response is any greater 
or any less in more general laminates, i.e., the (23/-60)«p 
laminate or designs 1 and 2, than in more conventional 
laminate, i.e., the symmetric quasi-isotropic design. It 
does appear that perhaps the magnitude of the response in 
the boundary layer can be controlled to some degree by 
laminate design, i.e., designs 1 and 2 vs. the symmetric 
laminate in fig. 31. The advantage, or disadvantage of the 
highly coupled nature of response of cylinders can only be 
assessed in the context of a specific application. 
However, it is clear that there is no limit to the type of 
response that can be achieved. 

Three directions for future research are clear. 
First, it would be useful to fabricate cylinders that 
exhibit a high degree of coupling and verify that the 
predictions are valid. The thermal response is felt to be 
particularly important because at the service temperature 
the cylinders are not at a constant radius toward the ends. 


This would be a problem in attaching the cylinder to other 
structural components. 

Second, the calculation of the response of a cylinder 
to nonaxisymmetric loadings is an important next step. 
Specifically, the response to bending loads is felt to be 
important enough to recommend it as a next step. 

Finally, determining the influence of geometric and 
perhaps material nonlinearities is considered important, 
particularly geometric nonlinearities. The radial 
deformation of the cylinder cannot exceed one-half to one 
wall thickness before geometric nonlinearities become an 
issue. Inclusion of these effects can be complicating but 
they are important at higher- load levels in both 
compression and torsion. Nonlinear effects in the thermal 
case are not felt to be as important. 



REFERENCES 


Hyer, M. W., "Some Observations on the Cured Shapes of 
Thin Unsymmetric Laminates," J. Composite Materials , 
vol. 15, Mar. 1981, pp. 175-194. 

Hyer, M. W. , "Calculations of the Room-Temperature 
Shapes of Unsymmetric Laminates," J. Composite 
Materials , vol. 15, July. 1981, pp. 296-310. 

Chaudhuri, Reaz A., Balaraman, K. , and Kunukkasseril , 
Vincent X. , "Arbitrarily Laminated, Anisotropic 
Cylindrical Shell Under Internal Pressure", AIAA 
Journal . Vol. 24, No. 11, 1851-1858, Nov. 1986. 
Sanders, J. Lyell, Jr., "Nonlinear Theories for Thin 
Shells," Quarterly of Applied Mathematics . Vol. 21, 

No. 1, pp. 21-36, 1963. 

Dym, Clive L. , and Shames, Irving H. , Solid Mechanics : 
A Variational Approach . McGraw-Hill Company, Inc., 
1973. 

Tsai, Stephen W. , Composites Design . Think Composites, 
1987. 

Donnell, L. H. , Beams. Plates, and Shells . McGraw-Hill 
Company, Inc., 1976. 

Reismann, H. and Pawlik, P. S., Elasticity. Theory and 
Application . John Wiley & Sons, Inc., 1980. 


